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INTEGRAL played a key role in discovering obscured sgHMXB in the Galaxy. We used XMM- 
Newton to perform X-ray wind tomography of a specific of these systems, IGR J 17252-36 16, fea- 
turing eclipses of the accreting pulsar. The X-ray band (0.2-10 keV) reveals vital information on 
the geometry of the surrounding gas probing simultaneously the absorption and the fluorescence 
emission. The XMM observations were scheduled to cover as many orbital phases as possible. 
Timing analysis allows the derivation of an accurate orbital solution and of the system parame- 
ters. Spectral analysis revealed remarkable variations of the absorbing column density along the 
orbit and of the Fe Ka fluorescence line around the eclipse. The combination of these observ- 
ables revealed a highly asymmetric and unprecedentedly extended structure in the stellar wind 
extending up to 2-3 stellar radii. The observations can be modeled in terms of three independent 
components: 

i) the unperturbed stellar wind 

ii) the contribution of a highly asymmetric hydrodynamic wind tail-like structure and 

iii) a cusp of material close to the neutron star. 

These dynamical structures are imaged for the first time in a sgHMXB and explain the source of 
the high obscuration. 
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1. Introduction 

High mass X-ray binaries (HMXB) consist of a neutron star or a black hole fueled by the 
accretion of the wind of an early-type stellar companion. Their X-ray emission, a measure of 
the accretion rate, shows a variety of transient to persistent patterns. Outbursts are observed on 
timescales from seconds to months and dynamical ranges varying by factors of 10 4 . The majority of 
the known HMXB are Be/X-ray binaries, with Be stellar companions. These systems are transient, 



featuring bright outbursts with typical durations on the order of several weeks [|16|]. A second class 
of HMXBs harbor OB supergiant companions (sgHMXBs) that feed the compact object by means 
of strong, radiatively driven stellar winds or Roche lobe overflow. Thanks to INTEGRAL, the 



number of known sgHMXB systems has tripled in the past few years [13]. 



IGR J17252 - 3616 was detected by ISGRI onboard INTEGRAL on February 9, 2004 among 



other hard X-ray sources []12[]. The source was first detected by EXOSAT (EXO 1722-3616) as 



a weak soft X-ray source, back in 1984 [15]. In 1987, Ginga performed a pointed observations 



and revealed a highly variable X-ray source, X1722-363, with a pulsation period of ~ 413.9 sec 



[JlOp. Additional Ginga observations revealed the orbital period of 9-10 days and a mass of the 
companion star of ~ 15 M e [^]. Both studies concluded that the system is a high mass X-ray 
binary (HMXB). 

INTEGRAL and XMM - Newton observations of IGR J17252 - 3616 allowed Zurita-Heras 
et al. (2006) to identify the infrared counterpart of the system, to accurately measure the absorbing 
column density, and refine the spin period of the system [fT7||. Thanks to the eclipses, an accurate 
orbital period could be derived from INTEGRAL data. Further RXTE observations helped identify 



a highly inclined system (i > 61°) with a companion star of M* < 20 M© and R* ~ 20 — 40 R Q [11 
Recent VLT observations help to infer the companion spectral type (between B0-B5 I and B0- 
Bl la) and radial velocity measurements [Bj. Its spectral energy distribution can be characterized 
by a temperature of T* ~ 30 kK and a reddening of Ay ~ 20 [Qj. 

2. Results 

The Science Analysis Software (XMM-SAS) version 9.0.0 1 was used to produce event lists for 
the EPIC-pn instrument. Spectra and lightcurves were built by collecting double and single events 
in the energy range 0.2 - 10 keV. The observations were scheduled to cover the orbital phases, 0.01, 
0.03, 0.08, 0.15, 0.27, 0.37, 0.40, 0.65, 0.79, and 0.91 inferred from the orbital solutions. 

Orbital solution 



We have obtained a refined orbital solution based on RXTE [11] and XMM-Newton Pulse 
Arrival Times [||] to derive accurate orbital phases. The orbital period is constrained using INTE- 
GRAL 2 long-term monitoring. We have obtained an orbital period of P = 9.742 ±0.001 days, a 
projected semi-major axis a x sini= 102 ± 8 lt-s, a mid-eclipse Tgo = M/D53761.67 ±0.1 and an 
90% upper limit e < 0.15 for eccentricity, yielding a mass function / ~ 11.7 M©. Corresponding 



'http://xmm.esac.esa.int/sas/ 

2 HEAVENS [114] interface available at www.isdc.unige.ch/heavens 
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masses of both donor star and compact object are in the range, Mob ~ 14 - 17 M© and M^s ~1.4 
-1.7 M©, respectively. 

Spectral Variability 

The spectra are always heavily absorbed below ~ 3 keV and contain an iron K-edge at 7.2 
keV. Throughout our observations an iron Ka line, at 6.40±0.03 keV, is always present. Some 
parameters (photon index, cut-off energy, blackbody templrature) did not vary significantly among 
the observations. To search for spectral variability we decided to fix them to their averages val- 
ues (Ec=8.2 keV, T = 0.02, kTgB=0.5 keV). The model was built by using an intrinsically ab- 
sorbed high energy cut-off power-law with a gaussian line, a black-body component responsible 
for the observed soft X-ray excess, and an overall absorption responsible for the interstellar ab- 
sorption fixed to the Galactic absorption N#=1.5xl0 cm" 2 . The model can be described as 
wabs* (bbody + gauss + vphabs*cutof f ) in XSPEC 3 using standard abundances ob- 
tained from [|l|]. All the normalization parameters and the intrinsic column density were free to 
vary. Figure 1 shows the variations of intrinsic N# (left) and the iron Ka line corrected equivalent 
width 4 (right) as a function of phase together with a wind model described in the next section. 
Huge variations of Nh can be seen, scaling up to a factor of 10. The unabsorbed 0.2-10 keV lumi- 
nosity level (outside the eclipse) is compatible with Lx ~ 10 36 erg s _1 , assuming a distance of 8 
kpc Minor variations on the flux can be interpreted as variation of the instantaneous accretion 
rate (M) onto the neutron star. Inside the eclipse the flux drops by a factor of ~ 200. In the follow- 
ing analysis, we assumed a distance D rs 8 kpc, a circular orbit (e = 0), and an edge-on geometry 
(i = 90°). 

3. Discussion 

Variability of absorption along the orbit 

We constructed ^ a 3D model of the OB supergiant stellar wind to investigated the behavior 
of the intrinsic column density, Nh, as a function of phase, and to identify potential structure of the 
wind during the orbit. We approximated the wind structure with two components, the unperturbed 
wind (Pwind) and a tail-like hydrodynamic perturbation (p ta a) related to the presence of the neutron 
star. These shocks are produced by hydrodynamical simulations [||] but produce a Nh of up to 
~10 22 cm~ 2 , which is too small to account for the variability observed in IGR J17252— 3616. 

The unperturbed stellar wind was modeled by assuming a standard wind profile [||] 



where V (r) is the wind velocity at distance r from the stellar center, is the terminal velocity of 
the wind, and j8 is a parameter describing the wind gradient. The conservation of mass provides the 
radial density distribution of the stellar wind. The unperturbed stellar wind is a good approximation 

3 heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 

4 As the X-ray continuum illuminating the gas emitting the Fe fluorescent line cannot be measured during the eclipse, 
we calculated a corrected Fe Ka EW by assuming a constant continuum flux of 1.8 x 10~ 3 ph keV -1 cm~ 2 s _1 
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Figure 1 : Left panel: Simulated N# profiles together with the data. Smooth stellar wind is illus- 
trated by cyan, green, and, red curves as a function of M/ 1)^=0.7, 1, 2x 10~ 16 M /km, respectively. 
Black solid line shows the total N# consisting of the unperturbed stellar wind (green curve) and the 
tail-like extended component. Right panel: Iron Ken corrected equivalent width during the orbit. 
Green curve indicated the perdition of the hydrodynamical tail and the black curve the addition of 
the central cocoon. 



within the orbit of the neutron star. Hydrodynamical simulations |Q] of HMXB have shown that 
the wind can be highly disrupted by the neutron star beyond the orbit. 

To estimate the terminal velocity of the unperturbed wind, we studied the Nh variability using 
three different sets of parameters (Fig. 1 left). The mass-loss rate and terminal velocity are con- 
strained by the data to be in the range M w / Voo ~ (0.7 — 2) x 10 M / km (j8 has a very limited 
impact on the results, so we used 0.7). 

The fraction of the wind captured by the neutron star could be estimated from the accretion 
radius r acc = 2GMx / (vf )rb + u 2 ) ~ 2 • 10 11 cm (where v or b = 250 km/s is the orbital velocity) as 
/ ~ xrl cc /4nRl rb ~ 7.5 • 10~ 4 . The mass-loss rate is therefore M w ~f- l M^\.5- 10~ 6 M Q /yr and 
the terminal velocity of the wind is constrained to be in the interval Voo ~ 250 - 600 km/h. 

In our simulation, we adopted a terminal velocity u»=400 km/sec, a stellar radius R if =29 R & , 
a wind gradient j3=0.7, and a mass loss rate M*= 1.35 x 10 -6 M Q yr _1 . 

We assumed that the tail-like structure is created very close to the neutron star and opens up 
with distance. The density of the material inside the 'tail' decreases with distance to ensure mass 
conservation. Its distribution follows a 'horn'-like shape with a circular section. We adjusted the 
density of the tail-like structure to match the observations. The density distribution p H .,-, 1( / + p ta u is 
displayed in figure 2 left panel. The supergiant is located at the center (black disk). The tail-like 
structure covers about half of the orbit. 

Figure 1 (left panel) displays the simulated Nh variability from the above density distribution 
together with the observed data points. The data and the model shows that the tail-like perturbation 
is essential to understand the observed variations. 
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Variability of the Fe Ka line during the orbit 

Assuming that the intrinsic X-ray flux is unaffected by the eclipse, the Fe Ka equivalent width 
drops by a factor ~ 10 during the eclipse in an orbital phase interval of ~ 0.1. This indicates that 
the radius of the region emitting Fe Ka is smaller than half of the stellar radius (< 10 cm) and 
far more compact than the tail structure responsible for the absorption variability profile. 

Outside the eclipse, the equivalent width of the Fe Ka line is of the order of 100 eV. Following 
Matt (2002) % and assuming a spherical transmission geometry, this corresponds to a column 
density of Nh ~ 2 • 10 23 cm~ 2 . As this additional absorption is not observed, the region emitting 
Fe Ka must be partially ionized. 

It is therefore very likely that the dense cocoon corresponds to the inner and ionized region 
of the hydrodynamical tail. We thus added this partially ionized cocoon in our simulations, using 
a density of 3 • 10 11 cm~ 3 within a radius of ~ 6 • 10 11 cm ~ 3R acc . The Fe Ka emissivity map 
(fig. 2 right panel) was calculated by applying an illuminating radiation field (~ 1/V 2 ) to the 
density distribution. Figure 1 right panel displays the resulting simulated profile of the Fe Ka 
equivalent width together with the observed data. The green curve shows the variations in the Fe 
Ka equivalent width expected from the wind density profile excluding the central cocoon, which 
obviously could not reproduce the data. The black curve accounts for the dense central cocoon. 
The exact profile of the eclipse is related to the size and density profile of the cocoon. No effort 
has been made to obtain an exact match to the data. 

The majority of the Fe Ka is formed in a region that is small enough to allow for pulsation 
of the Fe Ka line. We searched for such pulsations in our longest and almost uninterrupted obser- 
vation. Folded lightcurves were built in the energy bands 6.2 - 6.7 keV, 2 - 6.1 keV, and 6.8 - 10 
keV and resulted in a pulse fraction of 49 ± 5%, 58 ± 3%, and 57 ± 4%, respectively. The ratio of 
the line flux to the continuum in the energy range 6.2 - 6.7 keV is ~0.25. Assuming that the line is 
not pulsed, we infer an Fe Ka pulse fraction of ~ 50%, which is in good agreement with the above 
measurement. Weak Fe Ka pulsation can be explained if the cocoon is isotropic. 




-4-2 2 4 -4 -2 2 4 

R. R. 



Figure 2: Left panel: Number density distribution in the plane of the orbit including a smooth 
stellar wind and a tail-like perturbation. The black disk at the center represents the supergiant 
companion. Right panel: Integrated iron Ka emissivity (relative units) centered on the neutron star 
at phase (j) = 0.5. The extended tail can be seen on the right. 
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4. Summary 

The absorbing column density and the Fe Ka emission line show remarkable variations. The 
column density, always above 10 23 cm~ 2 , increases towards 10 24 cm~ 2 close to the eclipse, as 
expected for a spherically symmetric wind. The wind velocity is unusually small close to u M = 400 
km/s. An additional excess of absorption of 2 • 10 23 cm~ 2 is observed for orbital phases > 0.3, 
which is found to represent a hydrodynamical tail trailing behind the neutron star. 

During the eclipse, the equivalent width of the Fe Ka line drops by a factor > 10 indicating 
that most of the line is emitted in a cocoon surrounding the pulsar, with a size of a few accretion 
radii. This cocoon is ionized and corresponds to the inner region of the hydrodynamical tail 

The parameters of the IGR J17252 — 3616 are very similar to these of Vela X-l, except for 
the smaller wind velocity. We argue that the persistently large absorption column density is related 
to the hydrodynamical tail, which has been strengthened by the low wind velocity. The tail is a 
persistent structure dissolving on a timescale comparable to the orbital period. 

Our interpretation can be tested using numerical hydrodynamical simulations and high reso- 
lution optical/infrared spectroscopy. If confirmed, half of the persistent sgHMXB may have stellar 
wind speeds several times lower than usually measured. 
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